Helical Andreev bound states and superconducting Klein tunneling 
in topological insulator Josephson junctions 
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Currently, much effort is being put into detecting unconventional p-wave superconductivity in 
Josephson junctions based on topological insulators (TIs). For that purpose we propose to use su- 
perconducting Klein tunneling, i.e. the reflectionless passage of Cooper pairs through a potential 
barrier in a gated ballistic junction. This phenomenon occurs due to the fact that the supercurrent 
is carried by helical Andreev bound states (ABSs) characterized by spin-momentum locking similar 
to the normal-state carriers. We derive the spectrum of the helical ABSs and the corresponding 
Josephson current for a junction made on the surface of a three-dimensional TI. The superconduct- 
ing Klein tunneling is predicted to yield a non-sinusoidal current-phase relation and an anomalous 
critical current Ic that does not vanish with increasing barrier strength. We also analyze the depen- 
dence of the IcRn product (where Rn is the normal-state junction resistance) on the microscopic 
parameters of the superconductor/TI interface, which leads to lower IcRn values than expected from 
previous models of the proximity-effect Josephson junctions. 

PACS numbers: 



I. INTRODUCTION 

Unconventional p-wave superconductivity and 
Majorana states predicted in topological insulator 
(TI) /superconductor (S) junctionsi have become a 
topic of vigorous research (see e.g. reviews^ri). In 
three-dimensional (3D) TIs the superconductivity can 
be induced in the proximity to a conventional supercon- 
ductor (e.g. Al, W or Nb) deposited on the surface of the 
TI material^"— The unconventional superconductivity 
originates from the broken spin rotation symmetry of 
the helical surface states, allovi^ing for mixed singlet 
s-wave and triplet p-wave pair correlations, whereby the 
p-wave component inherits the spin-momentum locking 
of the surface states. Models for the superconducting 
proximit y effect in the TIs have been proposed in Refs. 
[l6l - [23[l4l l7l. including the case of disordered TIs in which 
the p-wave superconductivity is suppressedi^ 

Due to the induced p-wave superconductivity the TI 
Josephson junctions are expected to support gapless An- 
dreev bound states (ABSs)^ which have been studied 
theoretically in a variety of situations (see Refs. [2J- 
^M^ZM-M)- Transport mediated by the gapless 
ABSs may therefore serve as evidence for the unconven- 
tional superconductivity in TIs. Among such transport 
phenomena proposed in literature on TIs are the 47r- 
periodic Josephson effec t ^"^i^^" — and its signatures in the 
current noise^l the magnetically controlled Josephson 
effect r^i^i^^ the phase-dependent electric conductance 
and noised as well as the anomalies in the current-phase 
relation and Fraunhofer pattern^ 

Based on our recent proposal,— in this paper we study 
another transport manifestation of the induced p-wave 
superconductivity in TIs, a phenomenon analogous to 



Klein tunneling in normal helical metals. It should be ob- 
servable in ballistic TI Josephson junctions with a gated 
weak link [see Fig. [1^], which we model by introducing 
a potential barrier between the S terminals. We show 
that the induced p-wave superconductivity leads to he- 
lical ABSs that exhibit spin-momentum locking and de- 
pend strongly on the momentum direction with respect 
to the barrier. For the normal incidence at the barrier 
with 9 = Q [see Fig. [TJa], the ABSs are protected by 
time reversal symmetry and are gapless. In contrast, 
for the oblique incidence with 9 ^ Q there is no pro- 
tection since the incident and reflected trajectories are 




TI 




FIG. 1: (a) Schematic of S/TI/S junction with gated-induced 
barrier in the middle (see also text), (b) Scattering geometry 
in junction plane and polar plot of barrier transparency Tg as 
function of incidence angle 6 [see Eq. H3ip ]. 



not connected by time reversal symmetry. In this case, 
the ABSs acquire a minigap oc -/I — Tg depending on the 
normal-state barrier transparency Tq [see Fig. [T}d]. As 
a consequence, for the normal propagation the supercur- 
rent carried by the gapless ABSs is just the same as in 
the absence of the barrier. In this sense, there is a di- 
rect analogy with the normal-state Klein tunneling. Al- 
though the net supercurrent contains the contributions of 
all possible propagation directions, we demonstrate that 
the superconducting Klein tunneling is still identifiable 
by two features: the non-sinusoidal current-phase rela- 
tion and the anomalous critical current Ic that does not 
vanish with increasing barrier strength. We also analyze 
the IcRn product (where i?„ is the normal-state junction 
resistance) as function of the barrier strength and the mi- 
croscopic parameters characterizing the S/TI interface. 

The paper is organized as follows. In Sec. |TT]we review 
the superconducting proximity effect in a 3D TI, empha- 
sizing its two characteristic features: helical Bogolubov 
quasiparticles and the mixed s- and p-wave pair correla- 
tions. In Sec. mil we derive the helical ABSs and discuss 
their anisotropic properties. Section llVI presents the re- 
sults for the current-phase relation, critical current and 
the IcRn product. Section IVl summarizes our results. 



II. 



INDUCED SUPERCONDUCTIVITY IN TIs 



Model 



and Gz (ctq is the unit matrix), and v and /x are the sur- 
face Fermi velocity and energy, respectively). The ma- 
trix Hp is the representation of the surface-state Hamil- 
tonian in the four-component Nambu basis consisting of 
the wave functions ^-|-p, ^4.p, ^|_p and *J_p, where the 
star denotes complex conjugation. The tunneling self- 
energy E(e) ([2]) has off-diagonal entries describing the 
induced singlet superconducting pairing with strength 
A(e) and phase (f). The diagonal entries of S(e) account 
for the shift ir(e) of quasiparticle energies due to the 
tunneling. Both A(e) and r(e) are energy dependent, 
as they are related to the condensate and quasiparti- 
cle (momentum- integrated) Green's functions of the S, 
f,{t) = AJy/i^^Af and g,{e) = e/^e^-A^, where 
Ag is the energy gap in the S. 

The characteristic energy scale of S(e) is To which is 
determined by the tunneling coupling strength t and the 
normal-state density of states, A''^, in the S [see Eq. (|4])]. 
At high energies |e| ^ A^,, the induced pairing is sup- 
pressed: A(e) ^ FoAg/lel, whereas r(e) « Fq, yield- 
ing the quasiparticle escape rate into the S. In the zero- 
energy limit |e| — > 0, F(e) vanishes, and only the off- 
diagonal entries with constant A(e) sa Fg remain. In this 
limit the effective Hamiltonian ^ coincides with that of 
Ref. |l|. In our calculations we will use the more general 
Hamiltonian ([1]) given explicitly by 



h;" 



vcr • p — /i — iF(e) A(e)icrj 

— A(e)icrjye~' 



■5*0 



vcT* • p-l-^- iF(e) 



(5) 



We begin by setting up a microscopic model for the 
proximity effect in an infinite planar contact between 
a singlet s-wave superconductor (S) and a 3DTI. As- 
suming tunneling coupling between the systems, we will 
follow the general idea of McMillan's approacb^^, fur- 
ther adapted to various low-dimensional systems (see e.g. 
Refs. 36 - 39.16.17.40.411 . In this approach, the supercon- 
ducting proximity is accounted for by a tunneling self- 
energy S(e), leading to an effective Hamiltonian of the 
3DTI in the form 



H^J^ = Hp 



m 



S(6), 

-iF(e)cro 
— A(e)ii7j,e 



Hp = 





-h*_ 



-icf) 



A{e)iaye"t' 



A(e)=^Fo/,(e)=^Fo 
F(e)=Fo5s(e) = ro 



A.. 



v/i^^A|' 



v^ 



A2 



Fn = nt^N, 



(1) 
(2) 

(3) 
(4) 



In Eq. ([Il /ip = vcr • p — /i is the Hamiltonian of the 
surface state in the absence of tunneling (p — [pa:,py,0] 
is the carrier momentum on the surface, {h/2)<j is the 
spin vector expressed in terms of Pauli matrices <Tx,<Jy 



B. Helical quasiparticles, induced gap and mixed s- 
and p-wave pair correlations 

The p-wave correlations predicted in Ref. 1 for the 3D 
TIs arise from an interplay of the singlet s-wave pairing 
(X iuy and carrier spin helicity cr ■ p, where p denotes the 
unit vector in the momentum direction. The precise clas- 
sification of the pairing symmetry of the induced order 
parameter as well as the character of the single-particle 
excitations can be obtained from the Green's function Gp 
defined by the equation 



[6-i/;^^]Gp = /, 



(6) 



with I being the 4x4 unit matrix. We will consider Gp 
near the Fermi level positioned in the conduction band 
at high energies 



/i> A(e). 



(7) 



Inverting the equation for Gp we find the matrix elements 
of the Green's function as follows 



G, 



Gii|p Gi2|p 

G'21|p G22IP 



(8) 



where the indices 1 and 2 refer to the Nambu space, and 
the corresponding entries are 2x2 matrices in spin space 



given bjiisi 



Gulp 



Goain — TT 



'22|p 



1 (ctq + o- ■ p) [g(e) + v{p - pf)] 

1 -i(7y{aQ + a- ■ p)iay [£{e) - v{p - pp) 



Gl2\ 

G21IP 
where 



2 £^e) - v^p - pf)^ - A^e) 

1 (cro + cr-p)iCTj,A(e)e** 



P 2£2(e)_„2(p„p^)2„A2(e)' 
1 -z(T^(go + cr-p)A(e)e-»'^ 
2£2(e)_^;2(p_p^)2_A2(e)' 



(9) 
10) 

(11) 
(12) 



f(e) =e + ir(e), p = p/pF, <Ja = 



1 
1 



The single-particle excitations in the superconducting TI 
are described by the Green's functions Gii|p and G22\p- 
Both of them involve the projector ^(cto + cr • p), indi- 
cating that the quasiparticles are -1-1 eigenstates of the 
helicity cr ■ p, just like in the normal state (cf. Ref. \43) 
The pole of Eq. (j9|) yields the equation for the quasipar- 
ticle spectrum: 



S^(e)-v'ip-pFf-A'ie) 
For p close to pp the spectrum is 



± 



^v^{p~PF)^ + el. 



(14) 



(15) 



It features an induced energy gap Eg which satisfies the 
equation £'^{sg) = A2(£g). The latter can be explicitly 
written as (cf. Ref. mjt) 



1 - 



A. 



= 7 



1- 



£0 

A 



(16) 
where we introduce the dimensionless parameter 7 as- 
sumed small throughout. We may therefore seek the so- 
lution to Eq. (|16p in the form of the expansion in powers 

of 7: 



-^ = 7CI + 7^C2 -t- 7^C3 -I- .. 



(17) 



The constants ci, 02,03... are obtained by equating the 
coefficients at 7, 72, 7-^... on the left- and right-hand sides 
of Eq. p^ . Up to the cubic terms we find ci = — C2 = 1 
and C3 — 3/2, which yields the induced gap 



A. 



7- 



7^ + 27^ 



1 



7+27 



(18) 



This equation shows, in particular, that £g is reduced as 
the gap in the superconductor, Ag , becomes smaller (see 
also Fig. [2]). In Sec. IIVI we will discuss a similar depen- 
dence of the critical current of the S/TI/S junctions. It is 
also worth to note that the energy Fq can be expressed in 




As/i; 



(13) FIG. 2: Induced gap Eg (|18|l versus superconducting gap A3. 



terms of experimentally accessible parameters of a S/TI 
interface>i^ 



To 



hgn hv 
2^ A^' 



(19) 



where g„ is the normal-state interface conductance per 
unit area and kp is the Fermi wave-number of the TI 
surface state determined by the surface carrier density. 
The energy scale Fg can be extracted from the tempera- 
ture dependence of the critical current in short proximity- 
effect junctions. "^^ 

The off-diagonal entry G2i|p (fT2|) (as well as Gi2|p) is 
the Green's function of the induced superconducting con- 
densate. The symmetry of the induced pair correlations 
is identified by the spin structure of Eq. (IT^ : 



(iffy 



p)A(e). 



(20) 



It is a mixture of the singlet s-wave component (first 
term) and a triplet p-wave component (second term, see 
e.g. Ref. |4J). Equation ([20]) agrees with the results of 
Ref. Il6| for a large Fermi energy. The origin of the mixed 
s- and p-wave superconducting correlations is the bro- 
ken spin-rotation symmetry due to the helicity of the 
surface states. The situation reminds to some extent 
the mixed singlet-triplet intrinsic superconductivity pre- 
dicted for systems without inversion symmetry4^i^ The 
induced p-wave component inherits the spin-momentum 
locking (T • p of the normal-state carriers. Since it is an 
odd function of the momentum direction p, the p-wave 
component is suppressed in dirty TIs^"^ in which the elas- 
tic mean-free path i is much smaller than the supercon- 
ducting coherence length ^. In this paper we consider the 
opposite case of a clean TI defined by the condition: 



»^, ^^hv/VA'^{e)~£He). 



(21) 



In the next sections we discuss the link between the mixed 
s- and p-wave pairing (j20p and Andreev bound states 
(ABSs) in a ballistic S/TI/S Josephson junction. 



III. HELICAL ANDREEV BOUND STATES 

Let us consider a junction between two superconduct- 
ing TI terminals (see Fig. [T]). The junction length 
defined as distance between the terminals is assumed 
much smaller than the superconducting coherence length 
^ [see Eq. ^^]. Concretely, for Nb/TI contacts with 
typical parameters Ag « 1 meV, Fq « 0.2 meV and 
hv « 250 -;- 350 meV nm (see e.g. Refll4l). the junction 
length should be smaller than 1.25 -;- 1.75 /im, which is 
easily realizable experimentally. In this limit the junction 
can be described by an equation 



[H^ffir) + T,a,Uix)]^{r) - e *(r), 



(22) 



for the Nambu spinor vE'(r) in position representation. 
In Eq. dH]) H^ff{r) is the effective Hamiltonian © 
converted into position representation, whereby we in- 
troduce different phases 



(x < 0) = 0, 



0, 



(x > 0) = 0, 



(23) 



in the left (L) and right (R) terminals. In addition to 
that we introduce a potential U{x) = U5{x) in the mid- 
dle of the junction (the Pauli matrix Tz acts in the Nambu 
space). In practice, the potential barrier can be induced 
by a top-gate electrode as sketched in Fig. [1] The precise 
form of the scattering potential is not of crucial impor- 
tance here because we will focus on the effects arising 
from Klein tunneling that takes place independently of 
the barrier shape. We note that the related earlier stud- 
ies of Ref. il u28.2& assumed a scattering-free weak link. 
The influence of the normal scattering on the ABSs and 
the resulting superconducting Klein tunneling was first 
pointed out in Ref. 7, 

The ABS solutions of Eq. (|22p occur at energies e 
smaller than the induced gap and decay exponentially 
for X — > ±oo. We seek such solutions as the sum of 
independent modes with different ky wave numbers: 



j.kyy 



*(r) = V^fe {x)-^= 



(24) 



where ^^ [x) are evanescent functions for both a; > 
and a; < 0. For a large Fermi energy ^ the spinors 
^ky {x) are readily obtained from Eq. ([22]) as follows 



-^k{x>G) 



1 



a 



1 

— e' 

e^'^-a 

a 



R 



^i? 



fj>^ixkfcose-x/i (25) 



r^< ^—ixkp cos6 — x/^ 



-^kAx<Q) 



1 

e" 
-e'^a 

1 



"l 



^> ixkp cos 9+x/^ fOR] 



/^< ^—ixkp cos 9+x/^ 



where Cj^'^ are constants, the length-scale ^ is defined 
in Eq. (j2ip and other parameters and notations are ex- 
plained below: 



aiiW 



A(6)e- 



f(6)±*VA2(6)-f2(,)' 



cos^ = J\ — {ky/kpY 



(27) 



(28) 



Here a^ ^ (e) are the Andreev reflection amplitudes at 
the ends of the junction, the angle 9 indicates the parti- 
cle propagation direction (spanning the range from — 7r/2 
to 7r/2 around the a;-axis) for given open channel with ky 
smaller than the Fermi wave number kp. The matching 
condition for spinors (|25p and (P5| is obtained by inte- 
grating Eq. (|22|) around a; = 0, which yields 



*fe„M) = (l + iZr,a,)*fe^(+0), Z = U/hv, (29) 

where the dimensionless parameter Z characterizes the 
barrier strength. Inserting (|25|) and ([26l) into (|29p we 
arrive at four linear algebraic equations for the four con- 
stants C^ , C^ , C^ , and C^ . Equating the determinant 
of the system to zero, we then obtain after some algebra 
the following equation for the ABSs: 



T,[a+(e)-a+(e)][a]^(6)-aZ(e)], 



(30) 



where T{6) is the normal-state junction transparency: 

cos2(6l) 



To 



l-sin2(0)/(l + Z2)' 



(31) 



Using Eq. (P7)) for the Andreev amplitudes, we obtain 
from ((30)) an equation for e4i 




1 - Ta sin^ 



(32) 



Its solutions depend not only on the phase difference 0, 
but also on the relative energy Fp/Ag characterizing the 
S/TI interface. The latter dependence was not accounted 
for in the previous calculations of the ABSs in S/TI/S 
junctions. 

The solutions of Eq. (|5^ come in pairs e^ = ±|e|. In 
order to find |e| it is convenient to recast Eq. (|32l) as 



follows: 

2 



' 1 - T« sin-" ^ 



(33) 
where 7 is the small parameter defined in Eq. (J16p . We 
again seek the solution to Eq. (P5)) in the form of the 
expansion in powers of 7: 



a) A, 




J^=7/l(0)+72/2(<^)+7'/3(</') + 



(34) 



The functions /i(0), /2(0), /s (</>)■•■ are obtained by 
equating the coefficients at 7,7^,7^... on the left- and 
right-hand sides of Eq. (pS)) . Up to the cubic terms we 
find 



/3 = 



-/2 = 


l-Te 


. 2f 

sm - 


l-Te 


■ 2 4' 


1/2 1 

+ 2 



1/2 



(35) 



1 - To sin^ ^ 



3/2 



(36) 



Within this approximation the ABSs are given by^ 

^ll/2 



i?±(<^)«±rof (1-7 + 7') 

,-|3/2\ 



1 - T« sin^ 



T 
2 



1 - To sin^ ^ 



(37) 



The ABSs appear exactly below the induced gap given 
by(Ill,i.e. |i?±(0)|<£g = |i?±(O)|. 



A. Angle-selective protection of ABSs 

Figure [3] shows the ABSs for two distinct cases: nor- 
mal incidence at the potential barrier with 9 — (panel 
a) and oblique incidence with 6^0 (panel b). The per- 
pendicular propagation is refiectionless, with the trans- 
parency Tq — 1 for any barrier strength [see Eq. (j3ip and 
Fig. [TJd], which is the familiar Klein tunneling. In this 
case, the ABSs are gapless (see also Fig. [5^) 



E^iq^)^±ro 



(I-7 + 7') 





COS- 


4 


COS- 


3" 



.(38) 



We recall that gapless ABSs also appear in conventional 
ballistic Josephson junctions.™ Unlike the latter, the so- 
lutions (|38|) are protected by time-reversal symmetry. In 
order to check this we notice that the ± states in Eq. 
([55| are equivalent to two 47r-periodic branches crossing 
in the middle of the gap: 



E, 



(1,2) 



(</>) 



iFn 



9n <^ 7^ -i 4) 

(1 - 7 + 7^) cos I -H -^ cos^ I 

n(l,2). 



,(39) 



where we introduced notation E^ ' ' (0) to distinguish the 
47r-periodic solutions (IMl) from ([55)) . Let us now examine 



b)A, 




-^^ 



FIG. 3: Andreev bound states (|37p versus phase difference 
<j> for (a) normal incidence at potential barrier with 6 — 
and (b) oblique incidence with 9 = 7r/4; 7 = 0.3, Z = 0.5. 
Arrows indicate the 47r-periodic branches (|39p . Dashed line 
corresponds to fixed energy (|40|l . 



their wave functions at the same fixed energy (shown by 
the dashed line in Fig. [3^), which implies 



(40) 



E^'\^)^El^\2n-<j^) 



The wave functions of these states [see e.g. Eq. ([25]) 
for 9 = 0] appear to consist only of one-way propagating 
states: 



^i'\x) 



¥^\x) = 



1 
1 

-a+(0) 
a+(0) 

1 

-1 
a^{2TT — 






^^—ihpx — xj^ 



(41) 



(42) 



where l/2v^ is the normalization factor. Note that 
for the given energy (j40p the Andreev amplitude a^ in 

^Q (x) has the phase 27r— 0. It is now easy to see that the 
states (|41|) and (j42|) are connected by the time-reversal 
operation: 



v&f(x)=rozcT,*W*(^). 



(43) 






2.5e-05 


/^ 


\ 


-TC 


-71/2 


/ 


\ 


\ 


/ 


7t/2 


71 


\ 


vy 


-2.5e-05 





FIG. 4: Current-phase relation I{(j}) (|46p for different barrier 
strengths Z. Inset: current-phase relation /'(</>) (|47p for a 
reference conventional S/N/S junction for strong barrier with 
Z = 100. In both cases 7 = 0.2 and H = O.OlFo. 



the ABSs with oblique incidence at the potential barrier 
since in this case the incident and reflected trajectories 
are not related by the time reversal (see Fig. [Sb). The 
reflection from the potential barrier generates a minigap 
ex -^/l — Tg at (J3 — TT, Stt, ..., similar to the ABSs in con- 
ventional Josephson junctions4S 

Although ioT 6 y^ the ABSs are not protected 
from potential scattering, they still feature the spin- 
momentum locking and, therefore, can be called helical 
in the same sense as the normal-state 2D surface car- 
riers. It is worth noting that despite the dependence 
on the normal-state transparency Tg the helical ABSs 
([37]) reflect the nature of the induced superconducting 
condensate [see Eq. (^0]) ]. as they mediate Cooper-pair 
transport with the mixed s- and p-wave correlations. 



IV. CURRENT-PHASE RELATION I {(f)), 
CRITICAL CURRENT I^ AND IcRjv PRODUCT 



where tq denotes a unit matrix in Nambu space. More- 



,(1) 



,(2) 



over, vI/q (x) and ^q (x) are the eigenstates of the he- 
licity matrix Tza^'- 



In the short junction case the Josephson current is car- 
ried mostly by the ABSs and can be calculated in equi- 
librium using the formulai^ 



r.a,*|,'''^(x) 



±*, 



(1.2) 



(x). 



(44) 



This shows that the spin of these states is tied to the 
momentum direction. As a consequence of Eqs. (j43p 
and dBl), the ABSs (jH]) and ([42]) are orthogonal to each 
other and, therefore, immune to spin-independent poten- 
tial scattering. 

Clearly, the topological protection does not hold for 



ir/2 

7(0) = — — / d0cos0 



dEJW ,. Ej{4>) 



(45) 



where N = kpW/T: is the number of the open channels 
and S is the temperature in energy units. For the ABSs 
in Eq. (|37l) the current (gH) is 



J 



n'^) 



7r/2 



smc 



dO COS 



Tgll-.^t^-t^TgS^ll^^ 



(1-7 + 7^) 



1 - Tg sin^ 



1/2 



+ b' 



1 - Tg sin^ I 



3/2 



1 - Tg sin^ 



1/2 



2s/ro 



(46) 

where Iq = eNTo/2h. It is instructive to also consider a reference S/N/S junction of the same geometry, where N is a 
conventional 2D system with a spin-degenerate parabolic dispersion. In this system the proximity effect is described 
by the same self-energy as in Eq. ([2]). For the Josephson current in the reference S/N/S system we introduce notation 
/'((/)). The latter is given by the equation 



J'(<A) _ sin0 r'[l - 7 + §7^ - |7^T'sin^ |] ^^^^ 



(1-7 + 7^) 



1 - T' sin^ 



1/2 



+ h' 



1 - T' sin^ ■ 



3/2 



Il> 



1 - T' sin^ 



1/2 



2s/ro 



(47) 



where Iq = 2/o includes the factor of two due to the spin degeneracy and T' — 1/(1 + Z^) is the normal-state 
transparency of the N region with the same (delta-like) potential in the middle. For the non-helical carriers in the N 
neither T' nor I'{4>) is angle dependent. 



The equilibrium current-phase relation (|46p is 27r- in Fig. |4l Intcrcstigly, /(0) remains non-sinusoidal even 
periodic. It is plotted for different barrier strengths Z for very large Z. In contrast, for the conventional current 




FIG. 5: Normalized critical current versus barrier strength Z 
for S/TI/S and conventional S/N/S junctions; 7 = 0.2 and 
H = COIFq. Inset: Normal-state conductance versus barrier 
strength Z for TI and N channels (see also text). 



(|T7)) for Z ^ 1 we get /'(0) oc sint/), shown for compar- 
ison in the inset of Fig. |4l This difference reflects the 
hehcal character of the ABSs (1571) most pronounced in 
the gapless case = 0. 

Next we analyze numerically the critical current Ic de- 
fined as the maximum of /(</)) (j46l) at fixed parameters Z 
and 7. Figure [5] shows the normalized /c as function of 
the barrier strength Z. Initially decreasing with Z, the 
critical current saturates at a constant value for Z > 3. 
This is in sharp constrast with the conventional S/N/S 
junction for which I'^{Z) is strongly suppressed for Z > 3. 

The large supercurrent despite the presence of the 
strong barrier is a manifestation of the helical ABSs ([57)1 
that are able to transport Cooper pairs without reflection 
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FIG. 7: IcRn product in units of TrFo/e versus normalized 
superconducting gap A^; H — O.OIFq. 



for 9 = or with weak scattering ioi 9 j^ across the 
junction. This can be viewed as a superconducting ana- 
log of the normal-state transport dominated by the Klein 
tunneling. For comparison, the inset of Fig. [5] shows the 
normal-state (Sharvin) conductance of the TI channel: 



tt/2 

Gn = Go I d9cos9Te 




Go 






(48) 



along with the conductance of the conventional N channel 
G; = G'oT', with G'o = 2Go. 

Figure [5] shows another experimentally relevant char- 
acteristic, the product of the critical current Ic and the 
normal-state resistance i?„ = 1/G„ ([48|) normalized by 
ttTq/tt. In accord with the above discussion, the IcRn 
product is larger for the S/TI/S junction compared with 
the conventional S/N/S system. 

Finally, in Fig. [7] we plot the dependence of the nor- 
malized IcRn product on the superconducting gap A^ . 
The limit A^/Fo — > 00 corresponds to a frequently used 
simplified model of the proximity effect in which the in- 
duced pairing potential is approximated by an energy- 
independent constant A = eg = Fq. In this limit and for 
small Z we recover the known result 



IcRn ~ gFo/tt 



(49) 



However, as shown in Fig. [71 for any finite Ag/Fo the 
IcRn product is smaller than (|49|) . Thus, the dependence 
on the microscopic parameters A^ and Fq may be one 
of the factors contributing to lower values of the IcRn 
product in proximity-effect Josephson junctions. 



SUMMARY 



FIG. 6: IcRn product in units of vrFo/e versus barrier strength 
Z for S/TI/S and conventional S/N/S junctions; 7 — 0.2 and 
H = O.Oiro. 



Based on a microscopic model for the superconduct- 
ing proximity effect, we studied the Josephson effect in 
a short ballistic junction made on the surface of a 3D 
TI. The induced superconductivity in the TI is the mix- 
ture of singlet s-wave and triplet p-wave pair correlations. 



8 



Furthermore, the presence of the p-wave correlations is 
encoded in hehcal Andreev bound states characterized by 
the spin-momentum locking. We showed that the super- 
current carried by helical Andreev bound states decreases 
in the presence of the potential barrier in the junction, 
but cannot be completely suppressed due to the Klein 
tunneling of Cooper pairs through the barrier. Such su- 
perconducting Klein tunneling reveals itself in the non- 
sinusoidal current-phase relation and in the saturation of 
the critical current at high barrier strengths. This sug- 
gests that the p-wave superconductivity should be de- 
tectable in gated TI Josephson junctions. For that the 
gate-induced potential must be higher than the Fermi 



energy of the surface state. In this limit we also do not 
expect the shape of the induced potential barrier to play 
an essential role. 
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